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Abstract 

The sw(2)-algebraic model interacting with an environment is investigated from 
a viewpoint of treating the dissipative system. By using the time-dependent varia- 
tional approach with a coherent state and with the help of the canonicity condition, 
the time-evolution of this quantum many-body system is described in terms of the 
canonical equations of motion in the classical mechanics. Then, it is shown that the 
su(l, l)-algebra plays an essential role to deal with this model. An exact solution with 
appropriate initial conditions is obtained by means of Jacobi's elliptic function. The 
implication to the dissipative process is discussed. 
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1. Introduction 



One of the recent interests in nuclear theory is to investigate so-called thermal effect in 
nucleus and the dissipative process of the collective motion based on the quantum many-body 
theory or the quantum field theory. To deal with the thermal effects in the quantum field 
theory, it is well known that there are some methods such as the imaginary-time formalism, 
the real-time formalism in thermo field dynamics and so on. Four of the present authors 
(Y.T., J. P., A.K. & M.Y.) have proposed a method to describe the thermal effects and dissi- 
pation in many-body systems. B'B As for the description of the thermal effects in quantum 
many-particle systems, one of the key points in our formalism is to use the mixed-mode 
coherent state, which we have defined in the series of our papers, S^fl'i) as a trial state in the 
time-dependent variational approach. Furthermore, we learned that the su(l, l)-algebraic 
structure which the Hamiltonian has in the system under consideration plays an essential 
role when the non-equilibrium time-evolution and/or the dissipative process are realized in 
the quantum many-particle systems. Usually, so-called phase space doubling is carried out 
in order that the Hamiltonian has the su(l, l)-algebraic structure. & However, in general, 
the degree of freedom introduced by the phase space doubling cannot be interpreted as the 
environment or heat bath. We have given a formalism to describe the thermal effects and 
dissipation in a system in which the relevant and the irrelevant degrees of freedom interact 
each other and the irrelevant degree of freedom can be regarded as the environment. & Fur- 



find the su(l, l)-algebraic behavior in the system governed by the su(2)-algebra. The recent 



In this paper, we apply the above-mentioned formulation to the nuclear sw(2)-model 
which interacts with the environment represented by a harmonic oscillator. Main purpose in 
this paper is to show the utility of the su(l, l)-algebra. As was mentioned previously, with 
the help of the su(l, l)-algebraic structure, we will find the possibility of the description of 
the dissipative process. We are restricted ourselves in this paper to the zero temperature 
system. Namely, we will adopt a coherent state as a trial state instead of the mixed-mode 
coherent state in the time-dependent variational principle. Further, the exact solution for 
the time-dependent variational equation of motion will be obtained with the help of the 
su(l, l)-algebra. It is shown that the solution is given in terms of the elliptic function. This 
solution is similar to that encountered in the investigation of nuclear collective motions in 
the classical su(2)-models. B However, the physical situation is different. 

This paper is organized as follows. In the next section, we introduce the nuclear su(2)- 
model interacting with a harmonic oscillator. The su(l, l)-algebraic structure is found by 




presents the powerful technique to 



review of the series of our work is seen in Ref.§) 
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means of the Schwinger boson representation for the sw(2)-algebra. The coherent state, which 
presents the classical counterpart of the original quantum many-body system, is introduced 
in §3. By the use of the canonicity conditions, the state is parameterized in terms of the 
canonical variables. Thus, it is shown that the equations of motion derived by the time- 
dependent variational principle are nothing but the canonical equations of motion. In §4, 
the equation of motion is solved in a certain case, and it is shown that the solution is 
expressed in terms of Jacobi's elliptic function. The implication to the dissipative process is 
discussed in §5. The last section is devoted to a concluding remarks. 

§2. Nuclear sii(2)-model and the construction of su(l, l)-algebra 
Let us start with the following Hamiltonian : 



Here, a* and a are boson creation and annihilation operators, respectively, and the set 



In the case uo = x — 7 — 0, this model Hamiltonian is known as that of the pairing model 
in which the interaction is active for particle-particle pair with coupled angular momentum 
being 0. In the case oo = G = 7 = 0, this corresponds to the Lipkin model which consists 
of two energy-levels with the same degeneracy. The level spacing is 2he and the interaction 
is active for particle-hole pair with coupled angular momentum being 0. These are the 
schematic nuclear models governed by the s-u(2)-algebra. We call these models the nuclear 
s-u(2)-models. In the case G = x = 0, this model is nothing but the Jaynes-Cummings 
model if the representation of the sw(2)-algebra is adopted as the spin 1/2 one. The Jaynes- 
Cummings model is well known in the field of quantum optics. In general, we will regard in 
this paper the above model given in Eq.(2-1) as a schematic nuclear model including both the 
relevant and the irrelevant degree of freedom which are represented by the sw(2)-generators 
and a-boson, respectively. Namely, this model will be regarded as a typical model that the 
many-body system represented by the nuclear sw(2)-model interacts with the environment 
represented by a-boson. Other possibility to interpret this model is as follows : This model 
in Eq.(2-1) can be interpreted as a schematic nuclear model including both the collective and 




(2-1) 



(So, S±) composes the su(2) algebra : 



[a, a*} = 1 , 
[S+,S-] = 2hS ( 



[a, S ,±\ = , 
[S ,S ± ] = ±hS ± . 



(2-2) 
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the intrinsic degree of freedom which are represented by a-boson and the sw(2)-generators, 
respectively. 

In this paper, we use the Schwinger boson representation for the su(2)-generators.& The 
new boson operators (b, b*) and (c, c*) are introduced and the su(2)-generators are expressed 
in terms of these two kinds of boson operators : 

So = — (c*c — b*b) , S + = hcb , S- = hb*c . (2-3) 

The Casimir operator T su (2) is written as 

An (2 ) = S 2 + ^(S + S_ + S_S + ) = S(S + h) , (24) 

S= ^( c *c + b*b) . 

Further, the following operators can be also defined in terms of the original boson oper- 
ators (a, a*) and the Schwinger boson {b,b*) as 

h 

T = -{bb* + a* a) , T + = Tib* a* , T_ = Tiab . (2-5) 

These operators compose the su(l, l)-algebra, and the Casimir operator 1^(1,1) is also defined 
as follows : 

[f+, TL] = -2hf , [f , f±] = ±hf ± , (2-6) 
r s u(i,i) = T 2 - \{f + f_ + f_f + ) = f(f - K) , (2-7) 

f = -(bb* -a* a) . 

It should be noted that the two algebras are not independent in this model. Actually, we 
can calculate the commutation relations : [S, T ] = [S, T] = [T, So] = 0. However, we obtain 
[S,f ± ] ^ and [f ,S±] ^ 0. 

The Hamiltonian (2-1) is then expressed by using the su(2)- and the su(l, l)-generators 
together with their Casimir operators. The Hamiltonian is recast into 

H = H + H 1 + H 2 , (2-8) 
H = hu( a*a + -]+ 2eS - GS+S_ 



+2e(S -fo-f + h)- G(f + f )(2S - f - f + h) , (2-9) 
Hx = 7v^(aS + + S_a*) = 7v^(T_c* + cf+) , (2-10) 
# 2 = -f(S 2 +S!) . (2-11) 
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In any case, we have the commutation relation, [H, S] = 0. Thus, S is a constant of motion. 
Besides S, we have conserved variables in some cases : (i) If x — 0, that is, the case of no 
particle-hole interaction, then H does not contain S±. Thus, we can derive [H, T] = 0, so 
that T is a constant of motion in addition to S. (ii) If 7 = 0, that is, there is no interaction 
between the relevant and the irrelevant motion, then H does not contain any linear term 
with respect to a and a*. Thus, we obtain [H, Tq — T] = because To — T = ha* a. Namely, 
we have another constant of motion To — T in addition to S. (iii) If % = 7 = 0, then 
[H,T] = [H,Tq] = 0. Thus the variables T and To are conserved independently. 

In any case, the operators S, T, To and So commute each other. However, So is expressed 
in terms of the other three operators as So = S — To — T + H. Thus, it is enough to introduce 
the simultaneous eigenstate, \s,t,to), for S, T and To with the eigenvalues hs, M and hto, 
respectively. Except for the normalization constant, the eigenstate is given as 

\s,t,t ) = (a*)*-^*)* *'- 1 ^) 28 -'-* *^) . (2-12) 

Here, a\0) = b\0) = c|0) = 0. The eigenvalue equations are written as 

S\s,t,t ) = hs\s,t,to) , T\s,t,to) = ht\s,t,to) , To\s,t,to) = hto\s,t,t ) , 
S \s,t,t ) = h(s-t-t + i)\s,t,to) . (2-13) 

As for the su(l, l)-algebra, the case t > 1/2 will only be treated from now on. In this case, 
the eigenstate \s,t,to) is expressed by the use of S + and T + as 

\s,t,t ) = Af(f + ) t °-\S + ) 2s - t - t " +1 (b*) 2s - t " +t \0) 

Then, we can show the following relations : 

S-(b*) 2s - to+t \0) = f_(6*) 2s -* (,+ '|0) = . (2-15) 

Therefore, S + and T + play a role of the raising operators from the state (6*) 2s_ * 0+t |0). 



(2-14) 



§3. Construction of coherent state and classical counterpart 

In the previous section, we have introduced the eigenstate with respect to the operators 
S, T, T and ^o- Then, S + and T + are regarded as the raising operators. Thus, it is allowed 
to construct the following coherent state : 

|c) = iV c exp (vU^/h ■ f+) • exp (UYW^/h ■ S+) ■ exp ^ ■ WU^b*^ |0) . (3-1) 
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Here, V, W and Y are complex variables, U = J 1 + \ V\ 2 and N c is a normalization factor. 
The state |c) is rewritten as 

|c) = iV c exp (^Jj- Yc * + VU~ l b*a* + |0) 

= |c c ) <g> |c a6 ) , (3-2) 

,I|y|2 I /2, 



|c c ) = e-«' y ' -exph/-y C *| |0) , 



\ Cab ) = U- l e-*\ w \ 2 ■ exp [vU-^a* + \j^ WU ~ lb *j l°> • 

By introducing the following new operators, 

a' = Ua- Vb* , b' = Ub - Va* - ^2/Ti W , d = c - ^2/Ti Y , (3-3) 
it is shown that the state |c) is vacuum for a', b' and d : 

a'\c) = b'\c) = d\c) = . (34) 
Inversely, a, b and c are expressed as 



a = Ua' + Vb'* + y/2/h W*V , b = Va'* + Ub' + ^2/h WU , c = d + sj2/h Y . (3-5) 

By the use of the above relations, we can easily calculate the expectation values of the 
su(2)- and the su(l, l)-generators and their Casimir operators. For example, 

T=(c\f\c) = \W\ 2 + ~ , 

S = (c\S\c) = [\W\ 2 + \V\ 2 + \W\ 2 + \Y\ 2 , 

T = ( C |:fo|c) = [ \W\ 2 + |1 (l + 2\V\ 2 ) . (3-6) 



It is here noted that the dynamical variables are (V, V*, W, W*, Y, Y*). As is similar to the 
time-dependent Hartree-Fock theory with canonical form, & it is convenient to parameterize 
the coherent state in terms of canonical variables. For this purpose, we impose the following 
canonicity conditions : 

h 

(clihdfolc) = T - - , (c\ihd To \c) = , 
(c|z?ic^|c) = T — — , (c\ihd T \c) = , 

(c\ihdip\c) = S , (c\ihds\c) = . (3-7) 
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The sets of variables (Tq,(J)q), (T,(f>) and {S,ip) are those of the canonical variables. The 
above canonicity conditions can be easily solved by the use of the following relation : 



(c\ihd z \c) = i 



(\W\ 2 + ft/2) ■ (V*d z V - Vd z V*) 



+(W*d z W - Wd z W*) + (Y*d z Y - Yd z Y*) 



(3-8) 



As a result, the original variables are expressed by the canonical variables as 

V -- 



ITn — T 



2T 



e -*</>o e -«V'/2 



n 



Y = J-(2S - T - T + ft) • e~^ /2 



(3-9) 



Thus, we can obtain the expectation values of various operators with respect to |c) in terms 
of the canonical variables : 



S \c) = S-T -T + h , 



S+\c) = yj2(2S-T Q -T + h)JT - 1 lT^±L . P -i(4*«>)/2 



2T 



T)|c> =T , 



f+\c) = ^J(T - T)(T + T) • e ^e^ /2 , 
S + S-\c) = (2S - T - T + h)(T + T) , 
a*a\c) = (T — T)/K , 

2 



ft 

2T 

n 



s 2 + \c) = { l- ^ ] (t • T){2s v;, v • //.w ' 



1 - — I (T + T) (2S - T - T + ft)e i(<fe+0) 



Vha*S-\c) = ^(Tq - T)(T + T)(2S — T — T + ft) , 
v^a^+lc) = ^(Tq - T)(T + T) (2S — T — T + ft) e~^° . 



(3-10) 



Thus, the expectation value of the Hamiltonian with respect to |c) is given by 

H = (c\H\c) = H + H 1 + H 2 

H = (c\H \c) = u(T -T + ft/2) + 2e(S - T - T + ft) 
-G(2 1 S-T -T + ft)(T + T) 
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f(T ,T; S) , 



H, = (clH^c) = 2 7v /(T - T)(T + T)(2S — T — T + h) cos0 o 

= g(T ,T; S) cos 0o , 
#2 = (c|# 2 |c> = - X (l - V(2T))(T + T)(2S — T — T + ft) cos(0 o + 0) 

= /i(T o ,T;,S)cos(0 o + 0) . (3-11) 

Here, ifo> Hi an d -^2 have been defined in Eqs.(2-9)~(2-ll). 

The equations of motion are derived from the time-dependent variational principle : 

5 f t2 dt(c\ihd t - H\c) = . (3-12) 
Jti 

With the help of the canonicity conditions (3-7), the derived equations of motion have the 
form of canonical equations of motion : 
OH 

^0 = -TT- = g(T ,T; S) sin (f) + h(T ,T; S) sin((f) + (f)) , 
O(po 
dH 

^ = dT 

df(T ,T;S) dg(T ,T;S) dh(T ,T;S) 
= dn + 9Tq cos O + ^ cos(0 o + 0) , (343) 

r)H 

T = -— = h(T ,T;S)sm(<j )0 + <i ) ) , 
• OH 

df(T ,T;S) dg(T ,T;S) dh(T ,T;S) 
= + Qf cos 0o + 7^ cos(0 o + 0) , (3-14) 

dH 



if, 



dS 

df(T ,T;S) dg(T ,T;S) dh(T ,T;S) 

+ ^ cos 00 + ^ cos(0 o + 0) . (3-15) 



OS OS ^ u OS 

Let us consider simple cases : (i) If x — 0, then h(T , T; S) = 0. In this case, from Eq.(3T4), 
the quantity T is conserved in addition to S. This fact is, of course, originated from the fact 
[H,f] = mentioned in §2. It is thus necessary to solve the time-evolution of T : 

T = g(T ,T; S) sin 0o , 
■ df(T ,T;S) dg(T ,T;S) 

0o = QTq + QTq cos 0o . (3-16) 

The expectation value of the total Hamiltonian, which is a constant of motion, is expressed 
as 

H = /(T , T; S) + g(T , T; S) cos O = E . (3-17) 
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From Eqs. (|3-17|) and (3-16), we obtain 



*> ^' + (^J^m)' = l . (348) 



\g(T ,T;S)J \ g(T ,T;S) 
Namely, 

f 2 = [g(T , T; S)] 2 - [E - /(T , T; S)] 2 . (3-19) 

(ii) If 7 = 0, then g(T ,T;S) = 0. In this case, from Eqs.(3-13) and (3-14), the quantity 
To — T is conserved in addition to S. It is thus necessary to solve the time-evolution of 
T — T. However it is enough to know the time-dependence for T because T — T = : 

T = h(T , T; S) sin(0 o + 4>) , 
■ df(T ,T;S) dh(T ,T;S) 

0o = ^7 + 7^7 cos(0o + 0) • (3-20) 

The expectation value of the total Hamiltonian is expressed as 

H = /(T , T; S) + h(T , T; S) cos(0 o + <P) = E . (3-21) 
From Eqs.(3-20) and ( |3-21| ), we obtain 

I % \\(e -hto,t,s)Y = i 



\h(T ,T;S)J V HT ,T;S!) 
Namely, 

f 2 = [h(T , T; S)] 2 - [E - /(T , T; S)} 2 . (3-23) 
Our next task is to solve the equation of motion (|3 • 1 9|) or (|3-23|) which have the same form. 



§4. Utility of elliptic function for a solution of equations of motion 

In the previous section, it has been shown that in the simple cases the equation of motion 
for Tq which has a form in Eq.( ^-19| ) or ( |3-23| ) should be solved. Hereafter, let us consider 



the case (i), that is, \ = 0, in which the particle- hole interaction is not active. In the region 
T > T > h/2 and the adequate model parameters, the right-hand side of Eq. (|3-19| ) is always 
non-negative. However, we regard the right-hand side of Eq. ( p-19|) as a function of T with 
arbitrary parameter set : 

f 2 = F(T ) = [g(T , T; S)] 2 - [E - /(T , T; S)] 2 

= -G 2 (T - r ) (T - n)(T - r 2 )(T - r 3 ) , (4-1) 
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where we used the knowledge that the function F is the function of fourth order for To and 
the coefficient of T 4 is negative : F(Tq) = —G 2 Tq + • • •. It should be noted here that, in 
case (ii), if the parameters satisfy an inequality x 2 (l ~~ ^/(2T)) 2 < G 2 , the same situation 
as the case (i) is realized. From the fact F(T — > ±00) < 0, it is concluded that the two real 
solutions for F(T ) = except for T = must exist because the physical region F > must 
exist. In addition to this fact, in the second line in Eq.(4T), we assume that there exist four 
real solutions for F(Tq) = which are denoted as tq, ti, t 2 and T3 with To < T\ < t 2 < T3. 
We can numerically check the existence of four real solutions for F(Tq) = 0. For example, 
if the parameters are adopted as hG/e = 1/2, uo/e = 1, ^y/h/e = 0.1 and Eo/e = and the 
variables are S/h = 5 and T/h = 1/2, the above-mentioned situation is realized. Further, 
in this parameter set, the physical region is given in r < T < T\. Thus, hereafter, the 
physical region is taken as r < T < T\. Then, T corresponding to the physical value is 
parameterized as follows : 

T = r - ( To - n jx 2 , -1 < x < 1 . (4-2) 

Thus, we can derive the equation for x from Eq.(4T) : 

G 2 

x 2 = —(to - t 2 )(t - r 3 )(l - x 2 )(l - hx 2 )(l - k 2 x 2 ) , (4-3) 
where k-y and k 2 are defined as 

kl = ?o^J± } k2 = Z^Zi ; < k 2 < h < 1 . (4-4) 

TO — T 2 T — T 3 



Further, if we set 



*(*) = 1 ^ V, \ ' i-l<™<!) (4-5) 
yjl-k 2 + k 2 m 2 (t) 



the variable m(t) satisfies the following differential equation : 



m 2 = Z 2 (l-m 2 )(l-k 2 m 2 ) , (4-6) 
Z = -(J , (4-7) 

71 



where c^o, K and k are defined as 

^0 = ^ (to ~ t 2 )(t y - t z ) , 

f 1 dx 
K = / = K(k) , 

^0 ^J(l-x 2 )(l-k 2 x 2 ) 

k = 



2 



(r °- ri)(r2 - T3) (0<*<1) (4-8) 



l-k 2 \ (t ~ t 2 )(t x - r 3 ) 
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Here, K and k are the complete elliptic integral of the first kind and the modulus, respec- 
tively. If we take the initial value mit = 0) = 0, that is, To(t = 0) = tq, we thus obtain the 
following solution for m(t) : 

m{t) = sn(Zt, k) , (4-9) 

where sn(Zt, k) is Jacobi's elliptic function : 

/• m (*) dm 
Zt = I . (4-10) 

'(1 — m 2 )(l — k 2 m 2 ) 



As a result, we obtain T Q as a function of time t : 

rp ^ (r -r 1 )(ro-r 3 )sn 2 (Zt,fc) 
(ri - t 3 ) + (r - ri)sn 2 (Zt, k) 

It should be noted that, if the physical region is given by r 2 < T < r 3 , the solution for T 
is also expressed in terms of Jacobi's elliptic function similar to Eq. (|4-ll| ). 



In the Fourier series representation for Jacobi's elliptic function, the function k 2 sn 2 (Zt, k) 
is expressed as 

E 7r 2 00 n 

kV(Zt, k) = 1 - — - _ ■ g sinHn7[K{k , )/Km ■ cos(W) , (4-12) 



where k! = \J\ — k 2 and E represents the complete elliptic integral of the second kind : 



From Eqs.( [4-liD and ( |4-12|) , it can be seen that the time-dependent variable T has a funda- 
mental period to given by 

to = — ■ (4-14) 

UJ 



§5. Discussion 

As was mentioned in §1, it is interesting to investigate the dissipative process which is 
governed by the su(l, l)-algebraic structure. In this model-Hamiltonian we are treating in 
this paper, the su(l, l)-algebraic structure is hidden and appears by means of the Schwinger 
boson representation for sw(2)-generators. In the previous section, the parameter Xi which 
is the strength of the particle- hole interaction, is taken as 0. It has been shown that in 
this case there exists a fundamental period to in the solution expressed by the elliptic func- 
tion. However, if this period is long compared with the typical time scale we consider, the 
dissipative situation may be realized. In this section, we will investigate this possibility. 
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The fundamental period to is given in Eq. ( |4-14j ). Therefore, if the angular frequency 



ujq, which has been given in Eq.(4-8), is infinitesimal or very small value, the period to 
substantially becomes infinity. Thus, the periodicity disappears. In order to realize this 
situation, it is necessary that the modulus k is close to 1 because the complete elliptic 
integral of the first kind, K, reveals the behavior of a logarithmic divergence with k — > 1. 
As a result, o>o is close to and to becomes infinity. If T\ ~ T2, then k ~ 1 is realized. In this 
case, Jacobi's elliptic function sn(Zt, k = 1) is expressed in terms of the hyperbolic tangent, 
namely, 

sn(Zt, 1) = tanh Zt . (5-1) 
Thus, the solution T is expressed by 

TWo- ^T 3 +1 • (5-2) 

to— ti tanh 2 Zt 

This is a monotonically increasing function from tq at t = to n at t — > 00. Thus, the 
periodicity is lost. 

In this special case, the energy of the degree of freedom represented by a-boson increases 
monotonically, that is, 

E a = (c\huj(a*a + 1/2) \c) = u (t ~ T+ (5-3) 

increases for time because To increases monotonically. Here, it is possible to introduce a 
"coordinate" x which is defined and expressed in terms of the canonical variables : 

x 

V 'Zmu v 

x = ( c \x\c) = \It£-.\I^( W * V + WV *) 

— T rh-rh^ 
yflSS »- - - - » ■ (54) 

It is seen that the amplitude of x increases monotonically. Thus, x reveals the behavior of 
the amplified oscillation. If the degrees of freedom represented by a-boson and by the su(2)- 
generators are regarded as the environment and the relevant degree of freedom, respectively, 
the process in which the energy of the relevant motion dissipates to the environment is 
described. This is nothing but the dissipative process. This process may be interpreted as 
the "cooling" process for the relevant motion, although we do not treat the system at finite 
temperature. Actually, if the state |c) is replaced into the mixed-mode coherent state, the 
thermal effect may be treated in this framework. 
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§6. Concluding remarks 

We have investigated the nuclear sw(2)-model interacting with the harmonic oscillator. 
By the use of the Schwinger boson representation for the sw(2)-algebra, the su(l, l)-algebraic 
behavior was come in sight. In this model we have considered in this paper, it is allowed to 
obtain the exact solution for the classical equation of motion. It has been shown that the 
solution is expressed in terms of Jacobi's elliptic function. The solution has a fundamental 
period. However, if the period is so long, it may be possible to realize the dissipative process. 
In this paper, we are restricted ourselves to deal with the zero temperature case. If we 
apply our formalism to the finite temperature system, it is necessary to replace the coherent 
state into the mixed-mode coherent state. We have already investigated the mixed-mode 
coherent state in the many-body system composed of three kinds of boson operators. 0) The 
investigation of the behavior in the nuclear sw(2)-model interacting with the environment at 
finite temperature is left CIS db future problem. 
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